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Spacetime Emergence in the Robertson- Walker Universe from a Matrix model 
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Using a novel, string theory- inspired formalism based on a Hamiltonian constraint, we obtain 
a conformal mechanical system for the spatially flat four-dimensional Robertson- Walker Universe. 
Depending on parameter choices, this system describes either a relativistic particle in the Robertson- 
Walker background, or metric fluctuations of the Robertson- Walker geometry. Moreover we derive 
a tree-level Al-theory matrix model in this time-dependent background. Imposing the Hamilto- 
nian constraint forces the spacetime geometry to be fuzzy near the big bang, while the classical 
Robertson- Walker geometry emerges as the Universe expands. From our approach we also derive 
the temperature of the Universe interpolating between the radiation and matter dominated eras. 

PACS numbers: 98.80.-k, 11.25.Yb 



Recent astronomical data show that the current Uni- 
verse is very close to the spatially flat Robertson- Walker 
(RW) geometry Q . The Universe has evolved from a 
big bang singularity, near which quantum effects are 
expected to have played an important role. While a 
complete quantum gravity description of the big bang 
is unavailable, effective matrix model descriptions of 
string/ -theory on time dependent backgrounds have 
lead to a number of insights 0, [l]- However, the 
main focus of these studies so far has been on time- 
dependent and supersymmetry-preserving orbifold or 
plane wave backgrounds, but not on the physically rele- 
vant supersymmetry-breaking RW geometry. One tech- 
nical obstacle is that the latter lacks a null isometry. 
Hence the conventional light-cone quantization is not ap- 
plicable and a new approach is required. We develop such 
an approach in the present Letter: the characteristic fea- 
ture of our formalism is the presence of a Hamiltonian 
constraint, i.e. a vanishing energy constraint. Instead of 
fixing the light-cone momentum, we consider a sector of 
fixed Hamiltonian density. In this way, for the first time 
it becomes possible, at least at tree level, to construct an 
M2-brane or Al-theory matrix model [H for the realistic 
RW geometry and demonstrate the emergence of classi- 
cal spacetime from an originally fuzzy geometry. 

In this Letter we begin with the analysis of the geodesic 
motion of a single D-particle in the RW Universe. In 
particular, we propose a conformal mechanics model in- 
variant under one-dimensional diffeomorphisms. For two 
different parameter regimes and gauge choices, this me- 
chanical system describes either the geodesic motion of 
a point particle in the spatially flat RW background, or 
homogeneous metric fluctuations around the background. 
More precisely, in each case we flnd a conserved quantity 
and show that any sector of the fixed value of the quan- 
tity is described by the conformal mechanics. This is 
reminiscent of the AdS/CFT correspondence Q, where 
matter and gravity dynamics are mapped to each other. 



Here however both the descriptions of matter and gravity 
descend from the same CFT action. Finally, we derive 
a matrix model from the action of the bosonic M2-brane 
in the RW background, giving a many-particle gene- 
ralization of the single-particle action, as in flat space [3] . 
We show that imposing the Hamiltonian constraint in 
the matrix model ensures the emergence of spacetime. 
Emergence here means that the Hermitian matrices in 
the matrix model, whose eigenvalues encode the positions 
of D-particles, become simultaneously diagonalizable far 
away from singularities, such that the particle positions 
can be simultaneously measured, and geometric quanti- 
ties become classical 0] . 



I. CONFORMAL MECHANICS 

In ZJ-dimcnsional spacetime, requiring both homogene- 
ity and isotropy of the D—l spatial dimensions, the met- 
ric is constrained to the RW geometry [H 

ds^ = -dt^ + a^{t) [dx2 + k(x • dx)7(l - Kx^)] , (1) 

where k is a constant, and a{t) is the only undetermined 
scale factor depending on the cosmic time t. 

The conformal mechanics ( "CFT" ), which we will show 
to be closely related to the spatially flat (k= 0) RW Uni- 
verse, is of the general form 



CFT 



dt 



C2 



(2) 
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Here ci, ci are constants, ^pify is the only dynamical 
variable and riif) is the inverse of an einbein. Both 
and r\ transform under one-dimensional diffcomorphism 
t ^ s(t) as ((^(t),ry(i)) ^ {ip{s),^{s)|h). 

Integrating out the auxiliary variable r\ reduces the 
action to / A^P\J ci'.p^'^ + C2- On the other hand, fix- 
ing the gauge symmetry with an arbitrary function of 
time, r\ = f}{t), the mechanical system ^ essentially 
corresponds to the known conformal mechanics Q ■ The 
gauge fixed action is invariant under the transformation 
5tf ~ fi{f(p — \f'p): where f{t) is given by a solution 



mm 



0. This third order differential equa- 
tion has three sohitions which form the symmetry alge- 
bra so(l,2) [13. 

The Hamiltonian corresponding to the gauge fixed La- 
grangian is, with the canonical momentum p^p = f][t)ip, 



HcFT 



1 



m) 



(3) 



The physical states must lie on the surface of vanishing 
energy TicFT = in the phase space, as implied by the 
gauge fixing of diffeomorphism invariance. (Note that 
throughout the Letter, '=' denotes gauge fixings or on- 
shell relations.) 



A. CFT for Single D-particle Dynamics 

Our starting point is a novel formalism for a generic 
mechanical system satisfying the following two condi- 
tions: (i) The Hamiltonian is given by the inverse of the 
Lagrangian 



(4) 



where m is a constant mass parameter. This always holds 
for a square-root relativistic particle Lagrangian of the 
form C = —m^\ — gij{t, x)x^x^ , after the gauge fixing 
to identify the worldline affine parameter with the cos- 
mic time. (ii) There exists a conserved quantity with 
on-shell value v, such that for the sector of fixed v the 
Lagrangian is completely fixed on-shell as a time- and 
I'-dependent function, 



Then the square-root free Lagrangian 



:= 



£2 



2e^(i) e^t) 



(5) 



(6) 



together with the Hamiltonian constraint equivalently de- 
scribes the sector of fixed v (isj. This can be shown by 
observing that all the canonical momenta of (O take the 
same on-shell values as those of £. Further, the Hamil- 
tonian corresponding to £^ reads from ([4|) and ([6]) 



= {C^ -el)/{2e,). 



(7) 



The Hamiltonian constraint Tij^ = then implies ([5]). 
Henceforth, as an example (see also Q), we turn to a rel- 
ativistic particle in four-dimensional RW background ((!]) 
in spherical coordinates 



d,s2 



-Af- + a^{t) d/?2 + r2(p) (d6'2 + sin^^d^^) 



where r(p) = vx^ is given by sin(^/7;p) /^/k, p or 
sinh(-\/— K/j) I \J—K depending on n being positive, zero 
or negative, respectively. After identifying the world- 
line affine parameter with cosmic time, the point particle 



or D-particle [16| Lagrangian in the RW background be- 
comes 



m W 1 - {t) I + r2 (p)^^ + (p) sin^ 



(8) 



The canonical momenta for p, 6 and (j) are 

Pp ^ M{t)p, pg ^ M{t)r^p)e, p^ ^ M{t)r^{p)&\T?6^, 

(9) 

where we set M{t) := —vr?a}(t)IL for compact expres- 
sion. The corresponding Hamiltonian 



n = 



^m^-+ a-^{t) [p2 + r-2(p) (ji^ + pi/ sin'o'^ 



satisfies In spite of the arbitrariness of a{t), the 

dynamics is integrable, as there exist three mutually 
Poisson-bracket commuting conserved quantities 



- P^l sin 



P0 

2/1 



pl + jVr\p) 



constant , 
{vinnf' . 



(10) 

(11) 
(12) 



The dimensionless constant j plays the role of a classical 
so(3) angular momentum. Introducing a time dependent 
mass as an on-shell value of M{t), 



m^{t) m.a{t)^Ja'^{t) + v"^ = M{t) , 



(13) 



the Hamiltonian and the Lagrangian assume the on-shell 
values, Ti. = —nn?/C = m^{t)/a'^{t). As for ([6]) we have 



£. 



(p'+r' ip)e' + r' (p) sin'e <}>^) + 



{vm) 



2m At) 
(14) 

Again, all the off-shell canonical momenta of match 
with the on-shell ones of jS]). Further, the corresponding 
Hamiltonian 



2mAt) 



J2 



Pp 



— {vmY 



(15) 



exhibits the same mutually commuting conserved quan- 
tities as (fT0| - (fT2|) . Thus the surface of the vanishing 
energy = in the phase space of the dynamical 
system (|15p describes precisely the relativistic particle in 
the RW background for a sector of fixed v. Further the 
subsector of fixed angular momentum is reached by set- 
ting J2 = j(j+l), which reduces to the conformal 
mechanics ([3|) . In this way, the conformal mechanics ([3]) 
with the choice Lp= p, f]= m^(t), c\ ~ ~j(i+l) < 0, 
C2 = {ym)"^ > describes the geodesic motion of a rela- 
tivistic particle with respect to cosmic time in the spatially 
flat RW Universe, with fixed conserved quantities v, j . 



B. CFT for Homogeneous Gravity 

Pioline and Waldron Q observed that for a solely time- 
dependent, generic D-dimensional metric in longitudinal 
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gauge 



-e^g-^dt^ + p2/(^-i)^.^.dxMxJ' , det 5 1 , (16) 



the Einstcin-Hilbert action with cosmological constant 
reduces to a mechanical system for a relativistic "ficti- 
tious" point particle, 



-J d^x,/^ {R-2A) = ^y"di ^ 



lfD-2\ .2 ^ C A 



D-l 



Here V is the D—1 dimensional spatial volume, hence- 
forth normalized to one, and C ;~ |e~^£>^tr(g~^gg~^^) 
which contains a non-linear tr-model metric for the coset 
SLp - 1)/S0(L> - 1) as tr(5-\^g-i^) = hab{0)re''. In 
terms of the momenta = Li;g/(4e).pa = ^\Q^e~^habd^ , 

we identify C = ^h''\0)paPb as the kinetic energy on the 
coset space and write the Hamiltonian 



2e 

UJ 



ujC 



w := 8 



D-2 
D-l 



(17) 



C is conserved Q and, in fact, positive semi-definite pTj . 

The original motivation of Q to consider the ho- 
mogeneous modes only was based on the observation 
in [lO[ that near cosmological singularities, inhomoge- 
neous modes generically decouple. Here we make the fur- 
ther observation that the spatially flat RW geometry ([1]) 
is a special case of (|16p as e = p = a^~^, cjij = 6ij, and 
hence p6p is the most general homogeneous and non- 
perturbative fluctuation of the RW metric. In the cosmic 
time gauge e = g, the conformal mechanics with the 
choice (p= g, 77= ^uja^~^ or r] = u)/{A:tf), Ci = loC > 0, 
C2 = ujA describes the homogeneous metric fluctuations 
of the spatially flat RW Universe with respect to cosmic 
time, with fixed value of kinetic energy on the coset space. 
In particular, near the big bang (a~ 0), the choice of 
small f/ ~ ma\/ + v"^ and negative ci describes the 
matter, whilst large 7)= juja^~^ and positive Ci describes 
the gravity. As an example for the map between gravity 
and matter, the metric element g = det gij in grav- 
ity is mapped to the radial coordinate p of the particle 
trajectory. 



II. MATRIX MODEL FOR RW UNIVERSE 

Now we turn to the description of many D-particles 
in the spatially flat RW Universe. The description 
is gen erically given by a Yang-Mills quantum mechan- 
ics [ill, i.e. by a matrix model generalization of a single 
particle action. In a flat background, the coupling of 
the Yang-Mills potential [A^A-'p in the matrix model 
can be freely scaled and therefore is irrelevant. How- 
ever, in the RW background the coupling coefflcient is 
time-dependent, and cannot be simply deduced from 
the one-particle action . We determine this time- 
dependent coupling by deriving the tree-level A^-theory 



matrix model from the bosonic M2-brane action in the 
RW background. The dynamics of an M2-brane with ten- 
sion T embedded in a L'-dimensional target spacetime is 
governed by the Nambu-Goto action 



^1 



M2 



T / d^i 



y-det(9axfa^x''5^^(a;)) . (18) 



The three-dimensional worldvolume of the M2-brane is 
parameterized by coordinates a = 0, 1, 2, and the em- 
bedding functions are given by a;^(^), = 0, 1, • • • , D~\. 
Adopting the cosmic gauge i = — as well as the 
longitudinal gauge dtX^daX^ g^i^u = 0, a = 1,2 (see [l2| 
for a detailed procedure), the remaining worldvolume dif- 
feomorphism is, at this stage, static, i.e. ^"-independent. 
The Nambu-Goto Lagrangian in the spatially flat RW 
background then reduces to 



M2 



-Ta^[t)^(\-a'^(t)x^)d(tig , (19) 



where the determinant det^ — deXidaX^ dfiXi\ a, 13 = 
1,2, is taken over spatial M2-brane coordinates only. 
Spatial indices i,j are contracted with the flat metric 
5ij such that x^:= x'^x^ 5ij. With the momenta pi = 
Ta'^Xi ■\/detfJ/(l — a^i^), the equation of motion reads 

(20) 



Pi = da 



TaV(l-a^i2)det5 G' 



-la 



'pXi 



and the longitudinal gauge condition becomes 

PidaX'- = XidaX^ = . 



(21) 



In terms of the Nambu bracket {x, yj^^ := e^'^daxdpy 
{e^^ = 1, e"^ = — e^"),the determinant can be expressed 
as dctg ^ ^ {x\x^}^^{x,,Xj}^^. Further ^ and ^ 
imply respectively 

dt(p^) +T^a^dtdctg 



0. 



0. (22) 



In the following we consider the sector of solution-space 
with flxcd on-shcU value of the Hamiltonian density 



(23) 



Since Q transforms as a scalar density, fixing Q finally 
breaks the remaining static diffcomorphisms down to the 
static area- preserving ones. This gauge-fixed sector is 
then equally described by a square-root free Lagrangian 



Cn := i (r^a^ x'±i - fl-^T^a^ det G + n) 



(24) 



since the canonical momenta as well as the equations of 
motion are identical to the on-shcU ones of ([T9|) . Simi- 
larly to the one-particle case in section|ll the Hamiltonian 
constraint for Cn matches with ([53]) as 

Hn = {p^ + T^a^ det G - a^fl^) / {2a^n) = . (25) 

Identifying dct^ = m'^/{T^a'^) and n = ■m^{t)/a'^{t), 
the M2-branc Lagrangians and ([M]) would respec- 
tively reduce to the point particle ones ([5]) and for 



4 



K = 0. However, ([22]) then would imply that — 2m?a? 
is conserved and the Hamiltonian constraint (|25p could 
not be met. Therefore, unlike in the flat background [3| , 
in an expanding Universe the M2-brane dynamics cannot 
be consistently truncated to a point particle dynamics. 

The matrix regularization of (jM)) prescribes to re- 
place the dynamical fields a:;*(i,^") by time-dependent 
N X N Hcrmitian matrices X'^{t), the Nambu bracket 
{x%a;-'}^^ by a matrix commutator [l^, and 

the worldvolume coordinates ^" by non-dynamical matri- 
ces ^" satisfying the non-commutative relation [^", = 
ie"^ X constant [l^l ■ The resulting M2-brane matrix 
model is of the general form 



Cm2 = Tr 



(26) 



where Cl{t) := r2(t,^") and I is a constant length. The 
covariant time derivative DtX^ = X* — i[ylo,X'] in- 
volves a non-dynamical gauge field Aq, such that the 
matrix model admits a U(A^) gauge symmetry. With 
the canonical momenta = ^a?{VtDtX^ + DtX^^), 
5Aq gives the Gauss constraint for the physical states 
[P^^Xi] ~ 0, which is consistent with ([22|) . Combining 
the result of section |T] and with the identifications 
Q,{t) = ■m^{t)/a^{t), fj = m^{t) (|13p . we finally determine 
the precise form of the matrix model for N D-particles 
in the flat RW background, for which the Hamiltonian is 



-Tr 



Pf 



2a^{t)l~^[X\X^Y~{vmf . (27) 



One crucial difference between the two matrix models 
for the M2-brane (Uni) and for D-particles (gT]) is the 
last non-dynamical potential term which gives rise to a 



different Hamiltonian constraint, reflecting the different 
dynamical behavior of M2-brane and D-particles. We 
interpret the non-dynamical term as the temperature T 
of the Universe Q , since from ([77)1 and TYdo = this 
term corresponds to the average particle energy. For 
the D-particle case, we have T = mv^ / {2a\/a? + i^) 
which nicely interpolates between the temperatures of 
the early radiation dominated era, T oc and the late 
matter dominated era, T cx a~^. 

Both in the Hamiltonian corresponding to ([H]) 
and in ([27|l . the constraint 7i= gives rise to 
an explicit realization of spacetime emergence: 
-Tt[X\ X^f = Tt{X\ X'}\ {X,, Xj]t > must de- 
crease on-shell as the scale factor a(t) increases [isj . 
Near the singularity the X* do not commute and thus 
cannot be simultaneously diagonalizcd. This leaves 
the particle positions fuzzy, and hence the spacetime 
geometry they probe. As the scale factor grows or the 
Universe cools down, the fuzziness disappears and the 
classical Robertson- Walker geometry emerges. In 0], a 
similar scenario was found in plane wave backgrounds 
from the time-dependent coupling of the Yang-Mills 
potential. Our matrix model is, however, subject to the 
additional Hamiltonian constraint. 

As a generalization of the one-particle case in sectionU 
the inhomogcneous fluctuations in the RW Universe 
are expected to be mapped to the matrix model (|26p . 
Quantum corrections to our scenario can be calculated 
in analogy to 0], in particular they may restrict the 
time dependence of ait). 
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